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We study the generalized mutual information of the ground state of different critical quantum 
chains. The generalized mutual information definition that we use is based on the well established 
concept of the Renyi divergence. We calculate this quantity numerically for several distinct quantum 
chains having either discrete Z{Q) symmetries (Q-state Potts model with Q = 2,3,4 and Z{Q) 
parafermionic models with Q = 5, 6, 7, 8 and also Ashkin-Teller model with different anisotropies) or 
the U{1) continuous symmetries (Klein-Gordon field theory, XXZ and spin-1 Fateev-Zamolodchikov 
quantum chains with different anisotropies). For the spin chains these calculations were done by 
expressing the ground-state wavefunctions in two special basis. Our results indicate some general 
behavior for particular ranges of values of the parameter n that dehnes In- For a system, with 
total size L and subsystem sizes I and L — i, the In has a logarithmic leading behavior given 
by ^ log(^ sin(^)) where the coefficient c„ is linearly dependent on the central charge c of the 
underlying conformal held theory (CFT) describing the system’s critical properties. 

PACS numbers: 11.25.Hf, 03.67.Bg, 89.70.Cf, 75.10.Pq 


I. INTRODUCTION 

The entanglement entropy, as a tool to detect and clas¬ 
sify quantum phase transitions, has been playing an im¬ 
portant role in the last fifteen years (see [I| and references 
therein). In one dimension, where most of the critical 
quantum chains are conformal invariant, the entangle¬ 
ment entropy provides a powerful tool to detect, as well 
to calculate, the central charge c of the underlying CFT. 
For example, for quantum chains, the ground-state en¬ 
tanglement entropy of a subsystem formed by contiguous 
i sites of an infinite system, with respect to the comple¬ 
mentary subsystem has the leading behavior S = | In £ 
if the system is critical or S' = |log^, when the sys¬ 
tem is noncritical with correlation length ^ [3]. Although 
there are plenty of proposals to measure this quantity in 
the lab the actual experiments were out of reach so 
far. Strictly speaking the central charge of quantum spin 
chains has never been measured experimentally. Recently 
other quantities, that are also dependent of the central 
charge has been proposed [ 1 , 0 ] • Among these proposals 
interesting measures that, from the numerical point of 
view, are also efficient in detecting the phase transitions 
as well as the universality class of critical behavior, are 
the Shannon and Renyi mutual informations [MU (see 
also the related works [l0-[3). The Renyi mutual in¬ 
formation (the exact definition will be given in the next 
section) has a parameter n that recovers the Shannon mu¬ 
tual information at the value n = 1. The results derived 
in id indicate that the Shannon and Renyi mutual 
informations of the ground state of quantum spin chains, 
when expressed in some special local basis, similarly as 
happens with the Shannon and Renyi entanglement en¬ 
tropy, show a logarithmic behavior with the subsystem’s 
size whose coefficient depends on the central charge. 


Recently additional new results concerning the Shan¬ 
non and Renyi mutual information in quantum systems 
were obtained, see [l9l - l^ . There are also studies of the 
mutual information in classical two dimensional spin sys¬ 
tems [l^, [23 - [^ . It is worth mentioning that the Shan¬ 
non and Renyi mutual informations studied in the above 
papers, as will be defined in the next section, are ba¬ 
sis dependent quantities. It is important to distinguish 
them from the more known basis independent quantity, 
namely, the von Neumann mutual information. For re¬ 
cent developments on the calculation of the von Neu¬ 
mann mutual information in thermal equilibrium and 
non-equilibrium systems see 

Most of the results regarding the Shannon and the 
Renyi mutual information, except for the case of har¬ 
monic chains, are based on numerical analysis, especially 
for systems with central charge not equal to one. One 
of the main problems in a possible analytical derivation 
comes from the presence of a discontinuity at n = 1 of 
the Renyi mutual information. This discontinuity pre¬ 
vents the use of the replica trick, which is normally a 
necessary step for the analytical derivation of the Shan¬ 
non mutual information. 

In this paper we will consider, for many different quan¬ 
tum chains, another version of the mutual information, 
which is also parametrized by a parameter n that re¬ 
duces at n = 1 to the Shannon mutual information. The 
motivation for our calculations is two fold. Firstly this 
definition is more appropriate from the point of view of 
a measure of shared information among parts of a sys¬ 
tem, since it has the expected properties. This will be 
discussed in the Appendix. 

Secondly, this quantity does not show any discontinu¬ 
ity at n = 1, so it might be a good starting point for 
the analytical calculation of the Shannon mutual infer- 
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mation with some sort of analytical continuation of the 
parameter n. From now on we will call this new quantity 
generalized mutual information. 

Having the above motivations in mind we firstly calcu¬ 
lated numerically (using exact diagonalization) the gen¬ 
eralized mutual information for several critical quantum 
spin chains. We considered models with Z{Q) symme¬ 
tries like the Q-state Potts modes for (5 = 2,3 and 4, the 
Z(4) Ashkin-Teller model and the Z{Q) parafermionic 
models with Q = 5 — 8. We then calculated the gener¬ 
alized mutual information for quantum critical harmonic 
chains (discrete version of Klein-Gordon field theory) and 
also for quantum spin chains with 17(1) symmetry like 
the XXZ and the spin-1 Fateev-Zamolodchikov quantum 
chains. 

The structure of the paper is as follows: in the next sec¬ 
tion we will present the essential definitions of the Shan¬ 
non and Renyi mutual informations as well as generalized 
mutual information. In section three we will present the 
numerical results of the generalized mutual information 
for many different critical quantum spin chains. Finally 
in the last section we present our conclusions. 


II. THE GENERALIZED MUTUAL 
INFORMATIONS: DEFINITIONS 

Consider the normalized ground state eigenfunction of 
a quantum spin chain Hamiltonian \'4>g) = 'Yhi 6^" 

pressed in a particular local basis |/) = |ii, ^ 2 , • ‘ ‘ )i where 
* 1 , * 2 , • ■ ■ the eigenvalues of some local operators de¬ 
fined on the lattice sites. The Renyi entropy is defined 
as 

SK{X) = ^\nY,p^, ( 1 ) 

1 — n 

where pi = |a/p is the probability of finding the system 
in the particular configuration given by |/). The limit 
n —>■ 1 gives us the Shannon entropy Sh = — ^jpj\n.pj. 
Since we are considering only local basis it is always pos¬ 
sible to decompose the configurations as a combination of 
the configurations inside and outside of the subregions as 
|/) = |/a7a)- Oil® can define the marginal probabilities 
as PI A = PIaIa and pi^ = PIaIa ■ 

In a previous paper we studied the naive definition 
of the Renyi mutual information: 

IniA, A) = Sh^{A) + SKiA) - SK{A u A). (2) 

From now on instead of using pi^i^ we will use just p/. 
The known results of the Renyi mutual informations of 
quantum critical chains are obtained by using the defini¬ 
tion For special basis, usually the ones where part of 
the Hamiltonian is diagonal (see [HI), the definition ([2]) 
for the Renyi mutual information gives us a logarithmic 
behavior with the subsystem size, for arbitrary values of 
n. However, as observed numerically for several quantum 
chains (see [10, [ill El ) i shows a discontinuity at n = 1, 


that forbids the use of large-n analysis to obtain the most 
interesting case where n = 1, namely the standard Shan¬ 
non mutual information. Although the definition ([2]) has 
its own uses it is not the one which normally has been 
considered in information sciences. For example for 
n A 1 is not necessarily a positive function, a property 
that we naturally expect to be hold for the mutual in¬ 
formations. In this paper we consider a definition that 
is common in information sciences [Slj. The generalized 
mutual information with the desired properties, as a mea¬ 
sure of shared information (see Appendix), is defined as 

M- 


iniA,A) 


-- In y 

n — 1 ^ 


Pi 


n-l n-l ’ 

Pia 


I Pi A 


(3) 


where and pj ^, as before, are the probabilities that 
the subsystems are independently in the configurations 
\Ia) and |/^) that forms the configuration |/) that occurs 
with probability pi. 

Hereafter L will represent the size of the whole system 
and i and L — t the sizes of the subsystems. With this 
new notation one can write /„(A, A) as /n(^, L — £). This 
definition of the generalized mutual information comes 
from the natural extension of the relative entropy to the 
Renyi case and measures the distance of the full distribu¬ 
tion from the product of two independent distributions. 
In the limit n —>■ I one easily recovers the Shannon mu¬ 
tual information — = Sh{i) + Sh{L — tj — Sh{L), 

where Sh = — 'YiiPiAipi is the standard Shannon en¬ 
tropy. One of the important properties of J„, that is not 
shared by , is its nondecreasing behavior as a function 
of n (see Appendix). Our calculations for a set of distinct 
quantum spin chains will be done numerically, since up 
to our knowledge an analytical method to consider these 
quantum chains is still missing. 


III. THE GENERALIZED MUTUAL 
INFORMATION IN QUANTUM GHAINS 


In this section we will numerically calculate the 
ground-state generalized mutual information of two se¬ 
ries of critical quantum spin chains with slightly different 
structure. In the first part we will calculate the general¬ 
ized mutual information for systems with discrete sym¬ 
metries such as the (5-state Potts models with (5 = 2,3 
and 4, the Ashkin-Teller model and the parafermionic 
Z((5)-quantum spin chain [s^ for the values of (5 = 5, 6, 7 
and 8. In the second part we will calculate the general¬ 
ized mutual information for systems with 17(1) symmetry 
such as the Klein-Gordon held theory, the XXZ model 
and the Fateev-Zamolodchikov model with different val¬ 
ues of their anisotropy parameters. 





A. The generalized mutual information in quantum 
chains with discrete symmetries 
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In this subsection we will study the generalized mutual 
information of the ground state of different critical spin 
chains with Z{Q) discrete symmetries. The results we 
present were obtained by expressing the ground-state 
wavefunction in two specific basis where the systems 
show some universal properties. 



FIG. 1: The generalized mutual information 7„(f, L — £) of 
the L = 28 sites periodic Ising quantum chain, as a function 
of ln(Lsin(-^))/4. The ground-state wavefunction is in the 
basis where the matrices Si are diagoiral (S' basis). 


1. The generalized mutual information of the quantum 
Q-state Potts model and the quantum Ashkin-Teller model 


Our results show that the Q-state Potts model and 
the Ashkin-Teller model share a similar behavior. For 
this reason we discuss them together. The critical te¬ 
state Potts model in a periodic lattice is defined by the 
Hamiltonian 


L Q-l 

= + R't). ( 4 ) 

i—1 k—1 

where Si and Ri are Q x Q matrices satisfying the fol¬ 
lowing Z{Q) algebra: [Ri,Rj] = = 0 

for i j and SjRj = e^~S' RjSj and Rf = Sf = 1. The 
model has its critical behavior governed by a CFT with 
central charge c = 1 - where = 2 cos(;j^). 

The Q = 2 Potts chain is just the standard Ising quan¬ 
tum chain. The Ashkin-Teller model has a Z{2) 0 Z{2) 



FIG. 2; The generalized mutual information In{i,L — £) of 
the L = 28 sites periodic Ising quantum chain, as a function 
of ln(Lsin(^))/4. The ground-state wavefunction is in the 
basis where the matrices Ri are diagonal {R basis). 


symmetry and a Hamiltonian given by: 


L 

H = -Y, {s,S!+,+SfS,+i+AS?Sl,+R,+R^+ARl) , 

( 5 ) 

where Si and Ri are the same matrices introduced in the 
Q = 4 Potts model. The model is critical and conformal 
invariant for —1 < A < 1 with the central charge c = 1. 
It is worth mentioning that at A = I we recover the 
(5 = 4 Potts model and at A = 0 the model is equivalent 
to two decoupled Ising models. 

In the paper [llj| we already showed that the Shannon 
and Renyi mutual informations, as defined in ([ 2 ]), are ba¬ 
sis dependent. In other words one can get quite distinct 
different finite-size scaling behaviors by considering dif¬ 
ferent basis. Surprisingly in some particular basis, that 
we called conformal basis, the results shows some uni¬ 
versality. For example, the results for the (5-state Potts 
model and for the Ashkin-Teller model in the basis where 
the matrices Ri or the matrices Si are diagonal are the 
same, and follow the asymptotic behavior 


r T TvP 

/„(7,L-7) = ^ln(-sin(^)) + ..., (6) 

4 TT T 

with 


r I, n = l 


( 7 ) 


We should mention that in [13 , based on numerical re¬ 
sults, it was claimed that for u = I the coefhcient ci 
might not be exactly equal to the central charge. As it 
was discussed in [iSini it is quite likely that /„ is not a 
continuous function around n = 1 and so any attempt to 
do the replica trick using this definition of Renyi mutual 
information will be useless. This makes the analytical 
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calculation a challenge. This is an additional reason to 
examine the behavior of /„, besides being the correct ex¬ 
tension, from the point of view of a measure of shared 
information. Having this in mind we calculated the /„ 
for Q = 2, 3 and Q = 4 Potts chains and for the Ashkin- 
Teller model in the H and the S basis. We found that in 
some regimes of variation of the parameter n one can fit 
the data nicely to 

~ c Tj irf 

7„(£,L-^) = f log(-sin(-)) + ..., (8) 

4 TT L 

being c„ a monotonically nondecreasing function of n, 
consistent with what we expect for the mutual informa- 
tit see 

th 



FIG. 3: The ratio Cn/c of the coefficient of the logarithm 
in the equation ([8| and the central charge c for the Q-state 
Potts model with Q = 2, 3 and 4, and for the Ashkin-Teller 
model (A-T) with different anisotropies A. The Ashkin-Teller 
model at the isotropic point (A = 1) is equivalent to the 4- 
state Potts model. The ground-state wavefunctions are in the 
basis where the Si matrices are diagonal. The lattice sizes of 
the models are shown and the coefficients Cn were estimated 
by using the subsystem sizes i = 3,5, Int[L/2]. 

Here we summarize the results for the Q-state Potts 
and Ashkin-Teller quantum chains: 

1. The results in general depend on the basis we 
choose to express the ground-state wavefunction. 

2. The generalized mutual information follows ([ 5 ]) in 
the S and R basis but with different coefficients for 
different basis. To illustrate the logarithmic behav¬ 
ior we show in Fig. 1 and Fig. 2 the mutual infor¬ 
mation In for the Ising model (Q = 2) with L = 28 
sites and ground-state eigenfunctions in the S and 
R basis, respectively. We see, from these figures, 
that for subsystem sizes 7 > 3 we have the loga¬ 
rithmic behavior given by ([ 5 ]) up to n ~ 8 in the 
S'-basis and n Ri 4 in the i?-basis. As we can see 
our results does not exclude the existence of some 


□ les 



FIG. 4: Same as Fig. 3, but with the ground-state wavefunc¬ 
tions of the quantum spin Hamiltonians expressed in the basis 
where the matrices Ri are diagonal. The lattice sizes of the 
models are shown in the figure, as well as the subsystems sizes 
I used to estimate Cn- 

3. The coefficient of the logarithm c„ in ([8]) is a con¬ 
tinuous monotonically non-decreasing function of n 
and it follows the following formula in the S basis: 

Cn = cf{n), with /(I) = 1, (9) 

where c is the central charge and /(n) seems to 
be a continuous universal function independent of 
the model, as we can see in Fig. 3. In the case of 
the Ashkin-Teller model the results start to deviate 
around n = 6 from the ones obtained for the Potts 
models. As we can see in Fig. 3, the deviation point 
is dependent on the anisotropy parameter A of the 
model. 

4. In the case of the R basis , as one can see in Fig. 4, 
equation Q is still valid for values of n up to ^4. 
However the function f{n) is distinct from the one 
obtained in the S basis. As shown in Fig. 4, up 
to n = 2 the form of the function f{n) seems to 
be also independent of the model. This figure also 
shows that the Ashkin-Teller model has stronger 
deviations in this basis, as compared with the re¬ 
sults obtained in the S basis. In order to better see 
the difference of the coefficients c„ in the S and R 
basis, we present in Fig. 5 the data of Figs. 3 and 
4 for the Q = 2, 3 and 4 state Potts models. 

5. The coefficient of the logarithm in the S basis al¬ 
ways goes to zero as n —0, differently from the R 
basis where it approaches to a non-trivial number. 
This simply means that probably in the continuum 
limit all the probabilities in the S basis are positive 
but in the R basis some of them are zero. For the 
definition of the n = 0 case see the Appendix. 
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n 

FIG. 5: The values of the ratios c„/c of Figs. 3 and 4 for 
the Q = 2, 3 and 4-state Potts model are shown in the same 
figure, for comparison. 


Our numerical results indicate that Cn is a continuous 
function of n around n = 1. This means that /„ should 
be a continuous function with respect to n and so it is a 
better candidate to be used in techniques exploring the 
analytical continuation of the value n, as happens for 
example in the replica trick. However, the appropriate 
technique that may be used is still unclear to us. 

It is important to mention that the results obtained 
for the ratio c„/c in this section (Fig. 3) and in the sub¬ 
sequent ones (Figs. 8, 9, 11 and 13) are based on the 
linear fit with the ln[L sin(£7r/L)] dependence. These fit¬ 
tings were done by choosing a set of subsystems sizes. In 
all the presented figures we only depict results where a 
small variation of the number of subsystem sizes gives us 
estimated values of c„ that differs a few percent. As an 
example we consider the fittings obtained from the data 
of Figs. 1 and 2 for the Ising model with L = 28 sites 
and ground-state eigenfunction in the S and R basis, re¬ 
spectively. This is shown in Fig. 6. As we can see, while 
for the S basis the fitting is reasonable up to n = 8 in 
the R basis we do not have reliable results for n > 4. 


2. The generalized mutual information in the parafermionic 
Z{Q)- quantum spin chains 

In this subsection we consider the generalized mutual 
information for some critical spin chains with discrete 
Z{Q) symmetry and central charge bigger than one. The 
quantum chains we consider are the parafermionic Z(Q)- 
^antum spin chains [3^ with Hamiltonian given by [SJ, 

el 

L Q-1 

H = - ^ ^(5^53.-'= +i??)/sin(^fc/Q), (10) 

i=i k=i 



FIG. 6: The values of Cn/c obtained from the data of Figs. 1 
and 2 for the Ising quantnm chain with L = 28 sites and 
eigenfunction expressed in S and R basis. 



ln(Lsin(Til/L))/4 


FIG. 7: The generalized mutual information In{£,L — £) of 
the L = 10 sites periodic Z(7)-parafermionic quantum chain, 
as a function of ln(Lsin(^))/4. The ground-state wavefunc- 
tion is in the basis where the Si matrices are diagonal (S 
basis). 


where again Si and Ri are the Q x Q matrices that ap¬ 
peared in ([4]). This model is critical and conformal in¬ 
variant with a central charge c = 2{Q — 1)/{Q + 2). For 
the case where Q = 2 and Q = 3 we recover the Ising and 
3-state Potts model, and for the case where Q = 4 we ob¬ 
tain the Ashkin-Teller model with the anisotropy value 
A = ^. We have considered the models with Q = 5, 6, 7 
and 8 and the ground-state wavefunctions expressed in 
the S 01 R basis. The results for the several values of 
Q are shown in Figs. 7, 8 and 9. To illustrate the loga¬ 
rithmic dependence with the subsystem size £ we show in 
Fig. 7 in{£,L — (■), as a function of ln(Lsin(7r£/L))/4 for 
the Z{7) parafermionic quantum chain with L = 10 sites, 
with the ground-state wavefunction expressed in the S 
basis. In Figs. 8 and 9 we show the ratio c„/c of the log¬ 
arithmic coefficient of ([5]) with the central charge c for 
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the Z((5)-parafermionic models with ground-state wave- 
function in the S and R basis, respectively. The max¬ 
imum lattice sizes we used for the Z((5)-parafermionic 
models are L = 12,11,10 and 9 for Q = 5,6,7 and 8, 
respectively. The results we obtained are very similar 
to the ones we already discussed in the previous case of 
the Q-state Potts models. All the five properties that we 
discussed in that subsection are equally valid also for the 
Z((5)-parafermionic models. By comparing the results of 
Figs. 8 and 9 with Figs. 3 and 4 we observe that the 
function f{n) in (I5|) are quite similar for the two set of 
models, at least for values of n up to ~ 6. Probably the 
matching of these curves is not perfect due to the small 
sys- 



n 

FIG. 8: The ratio Cn/c of the coefficient of the logarithm 
in equation ([8l and the central charge c for the Z{Q)- 
parafermionic models with Q = 5,6, 7 and 8. The ground- 
states are in the basis where the Si matrices are diagonal. 
The lattice sizes of the models are shown in the figure and 
the coefficients Cn were estimated by using the subsystem sizes 
^ = 3,5, ...,Int[L/2]. 


B. The generalized mutual information of quantum 
chains with continuous symmetries 


In this section we consider the generalized mutual in¬ 
formation of critical chains having a continuous 17(1) 
symmetry. We studied a set of coupled harmonic oscilla¬ 
tors which gives a discrete version of Klein-Gordon field 
theory as well as the spin-1/2 XXZ and the spin-1 Fateev- 
Zamolodchikov quantum chains. The last two models are 
interesting since, like the Ashkin-Teller model, they have 
an anisotropy that gives us a critical line of continuously 
varying critical exponents but with a fixed central charge. 



FIG. 9: Same as Fig. 8, but with the ground-state wavefunc- 
tion of the quantum spin Hamiltonians expressed in the basis 
where the matrices Ri are diagonal. The lattice sizes of the 
models, as well as the subsystems sizes £ used to estimate c„ 
are shown. 


1. The generalized mutual information in quantum 
harmonic chains 

In this subsection we will first consider the generalized 
mutual information of the ground state of a system of 
generic coupled harmonic oscillators. Then at the very 
end we will confine ourselves to the simple case where we 
have only the nonzero couplings at the next-nearest sites, 
that in the continuum limit gives us the Klein-Gordon 
field theory. 

Gonsider the Hamiltonian of L-coupled harmonic oscil¬ 
lators, with coordinates t/i,..., and conjugated mo¬ 
menta TTl, . . . , TTi: 


n—1 


L 

E ' 

n,n' = l 


^n^nn' 4^n' 


( 11 ) 


The ground state of the above Hamiltonian has the fol¬ 
lowing form 

^ (deG^)i^_i<^l^i/2|^>^ (12) 

TT-^ 

For the general Hamiltonian (ED. one can calculate the 
two point correlators Xa = tr {pA4>i4>j) and Pa = 
tr (pATTiT^j) using the K matrix defined in dill) . The 
squared root of this matrix, as well as its inverse, can be 
split up up into coordinates of the subsystems A (size £) 
and A (size L — £), i- e.. 


K -^/2 


Xa Xaa\ 

.Ka Xa) 


1^1/2 




Here we chose the couplings so that we always keep the 
equalities = Paa and = Xaa- The spectra 
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of the matrix 2C = \/XaPa, can be used to calculate 
the Renyi entanglement entropy (see [ 13 ] and references 
therein) as 


Sn{i,L-i) 






In this formulation we only need the correlators inside 
the region A. Note that the above quantity is basis in¬ 
dependent and is considered as an usual measure of the 
quantum entanglement. Here we need to introduce this 
quantity just for later use. To calculate the generalized 
mutual information for a system of coupled harmonic os¬ 
cillators one first needs to fix the basis. Here we work 
in the position coordinate basis, however all the results 
are valid also in the momentum basis. One should notice 
that the same is not true if one works in a generic basis 
obtained through canonical transformations from the po¬ 
sition or momentum basis. In order to calculate In first 
we find p{^a) and as 


p(^’a) = 
p{^a) = 






(13) 

(14) 


where Pa = Pa - Paa{Pa) ^PaA ^nd Pa = Pa - 
PL{Pa)-^Paa- Since cj) takes continuum values one 
needs to consider the integral version of the equation ([3]) 
as follows: 


In = 


n — 1 


■ In / 


p"($) 


pn l($^)pri 1($^)’ 


(15) 


where p($) = iTop. Plugging Eqs. (IT^ . (fT3|) (fT4)) in the 
equation (1151) and performing the Gaussian integral one 
can derive the generalized mutual information 


where 


( XaPa XaaPa\ 
{x^aPa XaPaJ 



XaaPa^ 


• ( 21 ) 


There is an important remark that we should mention: 
in principle Eq. & makes sense only if n -|- (1 — n)T is 
a symmetric positive definite matrix. If we start with a 
symmetric positive definite matrix this is already 

warrantied for 0 < n < 1 but for n > 1 one needs to 
check its validity. This will be an important point when 
we study the short-range coupled harmonic oscillators. 
Finally one can write 


inii, L-i) = S2[P L-i) + Mn{l, L-e) = S 2 {i,L- i) 

- ^(^^_l) d^det(l - (1 - nfx\APAXAAPA) 


where XIn{i,L — £) is the only n dependent part. We 
notice here that by changing n to 2 — n we just change 
the sign of the second term, i. e., — P) = 

-Mn{PL-t). 



In = ^\n 


2(n- 1) 


det K 2 


det Pa det Pa 


In 


det — (n — 1) 

det 


The following determinant formulas 

det(PA) detP4 = detK^^^, (16) 

det(P^) detPA = detK^/^, (17) 

det Pa det = det Xa, (18) 

det Pa det = det (19) 

allow us to write 

in{PL-i) = S2{e,L-e) 


- r,/ ^ ' Indet(n-b (1 - n)T), (20) 

2(n — 1) 


FIG. 10: The second term in the equation (1221) . Mn{£,L — £), 
as a function of ln(Z/sin(^)) for periodic quantum harmonic 
chain with L = 120 sites. 

When n ^ 1 the second term vanishes and we recover 
the result of @ 


h{i,L-i) = S2{£,L-i). (22) 

For massless Klein-Gordon theory the above result in one 
dimension gives, as a consequence the well known result 
for the Renyi entanglement entropy [H, [131 , 

1 T ~p 

/i(£,L-£) = -ln(-sin(^)) + ..., (23) 

4 TT L 

where the dots are the subleading terms. Our numerical 
analyses indicate that for short-range quantum harmonic 
oscillators the matrix n -I- (1 — n) T is symmetric positive 
definite up to just n = ric = 2[42j|. The numerical results 
show that for the values 0 < n < 2 the equation ([5|) is a 
























very good approximation, as we can see for example in 
Fig. 10. The coefficient c„ of the logarithmic term in (|S]) 
is obtained from the fitting of the model with L = 120 
sites is shown in Fig. 11 and in the range 0.4 < n < 1.6 
surprisingly it follows the simple formula : 

Ti — 1 

Cn = f{n) = 1 + 4 , 0.4 < n < 1.6 . (24) 

This is the red line in Fig. 11. At n = 0 we expect 
zero mutual information for our system, this means that 
based on the symmetry n —>■ 2—n the coefficient for n = 2 
should be C 2 = 2. Finally one can conclude that for in¬ 
teger values of n = 0,1, 2 the coefficient of the logarithm 
is 

Cn = f{n)=n, n = 0,1,2. (25) 



FIG. 11: The coefficient of the logarithm c„ in the equation 
®. The lattice size L = 120 and the coefficients c„ were 
estimated by using the subsystem sizes i = 3,5,..., Int[I//2]. 
The red line is given by Eq. 


2. The generalized mutual information of quantum spin 
chains with continuous symmetries 

The Hamiltonian of the XXZ chain is defined as 

L 

Hxxz = <^J +1 + o-j o-J+i + Acr|cr|+i), (26) 

where cr^, cr^ and are spin-i Pauli matrices and A 
is an anisotropy. The model is critical and conformal 
invariant for —1 < A < 1 with a constant central charge 
c = 1, giving us a good example to test the universality of 
our results with respect to the change of the anisotropy. 
The long-distance critical fluctuations are ruled by a CFT 
with central charge c = 1 described by a compactified 
boson whose action is given by 

S =-^ J (j) = (j) + 2TTR, (27) 


where the compactification radius depends upon the val¬ 
ues of A, namely: 

R = yj— arccos A. (28) 

As it is shown in Fig. 12, in the basis, the generalized 
mutual information /„(£, L — i) shows the logarithmic be¬ 
havior given in ([81) only for n < 2. This can be simply 
understood based on what we observed for the chain of 
harmonic oscillators. One can look to the Klein-Gordon 
field theory as a non-compactified version of the action 
(EH). Since we showed that in that case the generalized 
mutual information is not defined beyond n = 2 we ex¬ 
pect the same behavior also in the compactified version. 
Note that in our numerical calculations one can actually 
derive spurious big numbers for the generalized mutual 
information even for n > 2, but we expect all of them 
go to infinity in the thermodynamic limit. This behavior 
seems to be independent of the anisotropy parameter A. 

The coefficient of the logarithm in (jS)) forn < 2 is again 
given by dH), as we can see in Fig. 13, with a function 
f{n) which fits to the results of the harmonic chain per¬ 
fectly. We also considered the results in the case where 
the ground state wavefunction is expressed in the cr^ ba¬ 
sis and, except around n = 1, the equation (I5|) is not a 
good approximation. The second 17(l)-symmetric model 



FIG. 12: The generalized mutual information L — i) of 
the periodic XXZ quantum chain with anisotropy A = —1/2, 
as a function of ln(sin( The ground-state wavefunction 

is in the basis where the af matrices are diagonal (ct^ basis). 
The results are for lattice sizes L = 28 and L = 30 and give 
an idea of the finite-size corrections. 

we considered is the spin-1 Fateev-Zamolodchikov quan¬ 
tum chain whose Hamiltonian is given by [38l| 

L 

Hfz = e^{a,-(a/)2-2(cos7-l)(a,^a/ + <a/-) 

+ (29) 

where S = are spin-1 SU{2) matrices, trf = 
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SiSi+i = cr^ + erf. The model is antiferro¬ 
magnetic for e = -fl and ferromagnetic for e = — 1. It 
has a line of critical points (0 < 7 < with a quite 
distinct behavior in the antiferromagnetic (e = -fl) and 
ferromagnetic (e = —1) cases. The antiferromagnetic 
version of the model is governed by a CFT with central 
charge c = | while the ferromagnetic one is ruled 
by a c = 1 CFT [i^. We calculated Inii, L — £) in both 
critical regimes where c = 1 and c = |, and for differ¬ 
ent values of the anisotropy. We found a very similar 
pattern as that of the XXZ quantum chain, as can be 
seen in Fig. 13. The equation ([5]) is valid for values of 
n < 2 and the coefficient of the logarithm follows ([9]) 
with a function /(n) which is quite similar to the one 
we found for the quantum harmonic oscillators and the 
XXZ Tn 

for c 



FIG. 13: The ratio Cn/c of the coefficient of the logarithm 
in equation ((S]) with the central charge c for the XXZ and 
for the spin-1 Fateev-Zamolodchikov quantum chains (F-Z). 
The XXZ (Fateev-Zamolodchikov) ground-state wavefunc- 
tion are in the {S^) basis. The results for the XXZ 
are for the anisotropies A = 0, —1/2 and in the case of 
the Fateev-Zamolodchikov model their are for the couplings 
7 = 7 r/ 3 , 7 r/ 4 . The lattice sizes of the models are shown and 
the coefficients Cn were estimated by using the subsystem sizes 
£ = 4,5,...,LI2. 


IV. CONCLUSIONS 

In this paper we calculated the generalized mutual 
information In{£,L — £), as defined in ([3]), for quan¬ 
tum chains describing the dynamics of quantum systems 
with continuous or discrete degrees of freedom. Most of 
our analysis was purely numerical due to the absence, 
at the moment, of suitable analytical methods to treat 
this problem. We considered several integrable quan¬ 
tum spin chains. These quantum chains either have 
a Z{Q) symmetry (like the Q-state Potts model with 
Q = 2, 3 and 4, the Ashkin-Teller model, and the Z{Q)- 


parafermionic model with Q = 5,6,7 and 8) or a 1/(1) 
symmetry (XXZ quantum chain and the spin-1 Fateev- 
Zamolodchikov model). We also considered the discrete 
version of the Klein-Gordon field theory given by a set 
of coupled harmonic oscillators. In this case we have 
a continuum Hilbert space. We observed that by ex¬ 
pressing the ground-state wavefunctions in general ba¬ 
sis the obtained results are distinct. However, similarly 
as happens for the quantity given in ([2]) (see [HI), 
our results on some special basis reveal some general fea¬ 
tures. These basis are the ones where the S' or i? oper¬ 
ators are diagonal, for the models with Z{Q) symmetry 
or the ones where or S^ are diagonal for the mod¬ 
els with t/(l) symmetry. In a continuum field theory 
description of these quantum chains these basis are ex¬ 
pected to be associated to the boundaries that do not 
destroy the conformal invariance of the bulk underly¬ 
ing Euclidean conformal field theory, and for this rea¬ 
son we call them conformal basis [ij. Our results indi¬ 
cate that in these special basis the mutual information 
has the same kind of leading behavior with the subsys¬ 
tem size £ as we have in the Renyi entanglement entropy, 
namely Ini£,L — /) ~ ^ ln(^ sin(^)), with a function 
Cn = c/(n), with /(I) = 1. Differently from the Renyi 
entanglement entropy where the equivalent function f{n) 
is universal (for any model and any basis) in the case of 
In our results indicate that the function /(n) depends on 
the special basis chosen to express the ground-state eigen¬ 
function of the particular model. For the set of Z{Q)- 
symmetric models we considered the function /(n), for 
n < 4, although different for the S and R basis are similar 
as the ones of the Q-state Potts chain {Q = 2,3,4) and 
the parafermionic Z{Q) quantum chains {Q = 5, 6,7, 8). 
In the case of the Ashkin-Teller model our results indi¬ 
cate that f{n), for n > 2, also depends on the anisotropy 
A of the model. On the other hand the models with 
continuum symmetry showed a similar behavior only for 
n < 2. For n > 2 we have strong evidences that most 
probably the generalized mutual information is not de¬ 
fined. It is quite interesting that in these cases one can 
understand most of the results by just studying simple 
short-range coupled harmonic oscillators. 

In order to conclude we should mention that an ana¬ 
lytical approach for the Shannon entropy or the Shannon 
mutual information (Ji or Ii in ([2) and (|3])) is a theoret¬ 
ical challenge. The analytical methods to treat this kind 
of problem normally use some sort of analytical continua¬ 
tion, in the parameter n, like the usual replica trick. The 
results we present showing the continuity of around 
n = 1, differently from what happens with J„, indicate 
that In is probably more appropriate for an analytical 
treatment. 
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V. APPENDIX: THE RELATIVE ENTROPY 
AND THE RENYI DIVERGENCE 

In this appendix we review the definitions of the rela¬ 
tive entropy and its generalization: the Renyi divergence. 
The relative entropy is defined as the expectation of the 
difference between the logarithm of the two distribution 
of probabilities p and q, from the point of view of the 
distribution p, i. e., 

D{p\\ q) = ^Piln—. (30) 

i 

It can be considered as a measure of the difference be¬ 
tween the two distributions p and q. Although it is not 
a symmetric quantity it helps us to define the mutual in¬ 
formation of the subsets X ans Y of the system as follows 

IiX,Y)=D{piX,Y)\\piX)p{Y)). (31) 

In words, the mutual information between two parts of 
a system is just the relative entropy between the distri¬ 
bution probability for the whole system and the product 
of the probability distributions of the different parts. It 
tells how much the different parts are correlated. The 
natural generalization of the relative entropy is the Renyi 
divergence and can be defined (see [3l| for example), as 

Du{p\\ q) = (32) 


It has the following properties: for n > 0 we have Dn{p || 
q) ^ Q and li p = q then we have Dn{p || q) = 0. The 
especial case n —>■ I gives the usual relative entropy. We 
also define the n = 0 case by: 


-Do(p II g) =-lng(i|pi > 0). (33) 


It is worth mentioning that using the above definition 
Dq{p II q) is not zero except when for all i’s for which 
qi > 0 also Pi > 0 holds. 

Another important property is the following (see [4l| 
and references therein): 

Theorem: the Renyi divergence is a continuous 

and nondecreasing function of the parameter n. 

Comparing (1321) with (l3T]) and (1301) the natural defini¬ 
tion of the generalized mutual information is 

/„(X, Y) = D„(p(A, Y) II p{X)p{Y)). (34) 

The above definition is different from as given 

by dl]), and has been frequently used in different areas of 
information science. 
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